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> 

^ ; Abstract 

q 

[~>^ . We analyse a unique set of time-resolved echelle spectra of the dwarf nova IP Peg, obtained 

q ; at ESO's NTT with EMMI. The dataset covers the wavelength range of 4000-7500i and 

shows Balmer, Hel, Hell and heavier elements in emission. IP Peg was observed one day 
after the peak of an outburst. 

The trailed spectra, spectrograms and Doppler maps show characteristics typical of IP 
rS ' Pegasi during the early stages of its outburst, such as prograde rotation from the accretion 

j^ ■ disc flow, chromospheric emission from the secondary and spiral arms. The high-ionisation 

line of Hell A4686 is the most centrally located line and has the greatest radial extension 
compared to the Hel lines. The Balmer lines extend from close to the white dwarf up to 
Ri 0.45/?L, with the outer radius gradually increasing when moving from H5 to Ha. The 
application, for the first time, of the modulation Doppler tomography technique, maps any 
harmonically varying components present in the system configuration, and this variability 
information is not considered in standard Doppler tomography studies. We find, as ex- 
pected, that part of the strong secondary star emission in Balmer and Hel lines is modulated 
predominantly with the cosine term, consistent with the emission originating from the irra- 
diated front side of the mass-donor star, facing the accreting white dwarf. For the Balmer 
lines the level of the modulation, compared to the average emission, decreases when mov- 
ing to higher series. Emission from the extended accretion disk appears to be only weakly 
modulated, with amplitudes of at most a few % of the non- varying disk emission. We find 
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no evidence of modulated emission in the spiral arms, which if present, is relatively weak 
at that our signal-to-noise ratio was good enough to put a lower detection limit of any mod- 
ulated emission at 5-6%. Only in one arm of the Hell A4686 line, is there a possibility of 
modulated emission, but again, we cannot be sure this is not caused by blending with the 
nearby Bowen complex of lines. 



1 Introduction 



Cataclysmic variable stars (CVs) consist of a Roche-lobe filling main sequence star, called 
the secondary star, orbiting a white dwarf - the primary star. Mass flows from the secondary 
star through the inner Lagrangian point towards the primary star, ultimately leading to the 
formation of an accretion disc (AD). After the AD's formation, a bright spot (BS) is formed 
at the contact point where the gas stream hits the AD. Dwarf novae (DNe) are a subclass of 
cataclysmic variables with two different brightness states, quiescence and outburst. In the 
latter, the system usually brightens by 2-5 mag with a time-scale that can vary from tens to 
hundreds of days. DNe can also show a third brightness state, the superoutburst, where the 
system brightens by typically an additional magnitude or two compared to normal outburst. 

IP Peg is the brightest eclips i ng (/ rj 81°) DN above the period gap, and was discov- 
ered by iLipovetsky & StepanianI (|l98lh . It has a 3.8-h orbital period and a V magnitude 
in the range of 14—12 mag, the extreme values corresponding to quiescence and outburst, 
respectively. 

From the first photometric light curves during quiescence and return to quiescence 



^rom tne tirst pnotometric lignt curves aurmg quiescence ana return to quiescence 
dWood & Crawford . 1 1 9861 : IWood et al.lll989l) . it was evident that the BS produces a promi- 



nent orbital hump, when it rotates into view from behind the disc with a single-stepped 
ingress and double-stepped egress superimposed on a more gradual disc-eclipse. Infrared 
light curves showed 0.2 mag ellipsoid al variations frorn the se c ondary star superpo sed on 
a deep eclipse of the WD and the BS dSzkody & Mated. Il986h . Iwood et all (|l989h . based 
on the analysis of 39 eclipses, found out t hat IP Peg has a r apidly increasing orbital pe- 
riod. However, a new study performed by IWolf et al.l (|l993h and combining 49 eclipses 
along with previously reported ones, reached the conclusion that the orbital period varies 
sinusoidally with a period of 4.7 years. This variation they attributed to a third body, a late 
M dwarf star with a most probable mass of 0. IOM0. They also showed that the rate of 
shrinkage in disc radius during decline from outburst is constant during diff erent outburst 
cycles, following an exponential decay. A study of the flickering in IP Peg (IBruchl . I2OO0I) 
showed that the BS is the dominant source of flickering, with the AD and WD contributing 
practically nothing. 



Spectroscopy of the red star (IMartin et all. 1 19891) showed TiO bands and the Nal dou- 
blet (red dwarf features), several broad double-peaked emission lines (disc features) and 
a nan^ow third peak in the emissi on lines moving in anti-phase with the double peaks (a 
chromospheric emission feature). iBeekman et al.l (I2OOO1) performed a TiO study of IP Peg 
and concluded that TiO instead of the Nal doublet should be used for radial velocity studies 
of the secondary star. In this way they found that the contamination from disc emission, 
which in other studies returned an elliptical orbit, could be avoided and higher signal-to- 
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noise ratio could be achieved. They gave a final projected orbital velocity of the compan- 
ion star K2 = 331.3kms^' and refined values of the primary and secondary star masses 
Ml = 1.O5M0 and M2 = 0.33 M©, respectively. Covering a complete cycle from quies- 
cence to outburst and quiescence again, the spectral characteristics of IP Peg will now be 
described. 

Quiescence spectra (IMarshl. Il988r) show broad double-peaked emission lines on a blue 
continuum, a strong orbital hump in the continuum (but much weaker in the lines) and a 
flat B aimer decrement indicative of optically thick emission. The trailed spectra show the 
signature of a prograde rotation of the AD, where the blue-shifted emission is eclipsed prior 
to the red-shifted one, a weak S-wave component in the red-shifted wing and a 10% blue/red 
asymmetry before eclipse, most pr obably arising from line emission from the BS. 

Early rise-to-outburst spectra dWolf et al.l 1 1998!) show Balmer emission, TiO absorp- 
tion bands, KI and Nal absorption dou blets, a weak S-w ave and two bright spots in both 
Doppler and eclipse maps. The study of IWolf et al.l (|l998l '). as well as a re-analysis of thei r 
results using an extension of the classical eclipse mapping method (IBobinger et al.l . 1 19991) . 
concluded that the first spot should be identified as the BS while two possibilities existed 
for the second spot. It either was the beginning of formation of a spiral arm or a sec- 
ond inner BS due to gas -stream overflow as predicted by, e.g., simulations performed by 
Armitage & Liviol (Il996l) . 



Rise-to-outburst spectra (IHarlaftis et al.l . ll994) show a decrease in the equivalent width 
of Ha emission line, down to a level similar to that observed immediately after outburst and, 
simultaneously, weak Ha emission from the secondary due to irradiation of the secondary 
star by the boundary layer. 

Near-outburst spectra (IPiche & Szkodyl. 119891) show unusually strong high-excitation 
lines, Balmer lines in emission, indications of a very disturbed AD with a non-uniform 
He emission distribution, an outflowing wind from the WD vicinity and chromospheric 
emission from the secondary. 



Outburst spectra as presented by Marsh & Homel (| 19901) show: i) strong Balmer and 
Hell A4686 double-peaked emission lines rising from the AD, ii) non-Keplerian emission 
filling in the centre of Hell A4686, probably due to a compact outflowing wind or inflow- 
ing magnetic accretion columns from close to the WD, iii) Balmer but not Hell emission 
from the secondary, located in the polar region on the side facing the AD, iv) no evidence 
for an enhanced gas-stream compared to quiescence, v) large azimuthal asymmetry in the 
outer disc, with the blue-shifted side of the disc being brighter and possibly more extended 
at phase and vi) a sharp transient com ponent in the He l l line possibly generated from 
the interaction of the gas stream and disc. iMarsh & Homd (Il990l) suggested that He emis- 
sion can be explained as photoionisation of the disc by the boundary layer and found that 
Balmer emission in outburst can be attributed by 3/4 to the disc, 1/4 to irra diation of the 



red sta r and also to boundary layer inadiation. Another study performed by ISteeghs et al. 



(119961) showed, from trailed spectra and Doppler maps, asymmetric disc emission, strong 
secondary star emission and a low velocity emission in the Ha line only. This stationary 
emission, appearing as an emission blob in the centre of the Doppler maps, was interpreted 
as slingshot prominences, material trapped in mag netic loops and co-rot ating with the sec- 
ondary star. During another outburst maximum, IHarlaftis et al.l (| 19991) found metal line 
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emission such as Mgll. Doppler tomography of the He and metal Unes locates those lines' 
emission on the inner Roche lobe of the s econdary s tar. 

Early decline-from-outburst spectra (IHessmanl. 1 19891) reveal for the first time in such 
strength a blend of high-excitation lines of Hell, CIII and NIII along with a very strong 
chromospheric emission-line comp onent from the second a ry. 



Decline-from-outburst spectra (IMorales-Rueda et al.l l2000f) show iiTadiation-induced 
emission from the companion star in Balmer, He, Mgll and CII lines. This emission, located 
only near the poles of the secondary, suggests shielding of the equatorial regions of the 
secondary star by the vertically extended AD. 



Early-quiescence spectra in the infrai^ed (|Littlefair et all l200lh reveal the existence of 
mirror eclipses in IP Peg. A mirror eclipse is an eclipse of the secondary star by an optically 
thin AD and appears in the trailed spectra as a reduction in the equivalent widths of the 
lines, at phase 0.5. A min^or eclipse also has a mirror-symmetry with the classical rotational 
disturbance, i.e. the red-shifted part of the line is eclipsed prior to the blue-shifted one. 

Even though IP Peg had shown many pecul iarities and noveltie s in the studies cited 
previously, it had to reveal another breakthrough. ISteeghs et all (| 19970 were the first to find 
convincing evidence for the existence of spiral structure in the AD of a close binary Spiral 
wav es had l ong b een predicted in studies of tidal interactions and numerical simulations 
(see iBoffinI (1200 lb and references therei n) and may play a c rucial role in our understand- 
ing of the angular momentum transport. ISteeghs et all (Il997h obtained spectrophotometric 
observations of IP Peg during the rise to an outburst, covering Ha and Hel A6678. Disc 
emission was centred on the WD and had a strong azimuthal asymmetry in the forms of 
two spiral arms. The repeat abihty of the two-arm ed spiral structure was confirmed during 
another outburst maximum. iHarlaftis et all (119991) applied Doppler tomography of the Hell 
A4686 line and clarified the shock nature of the spiral arms which showed a jump of a factor 
of more than 2 in intensity. Additionally they calculated an azimuthal extent of 90°, a 30° 
upper limit in opening angle and a 15% contribution of the shocks to the total disc emission. 
Morales-Rueda et al.l (120001) showed that the spiral pattern persists even 5 and 6 days after 
the outburst maximum, without showing any diminution in stre ngth. Time-resolved optical 
spectra obtained during late outburst by iBaptista et al.l (120001) reveal, through the eclipse 
mapping technique, that spirals are present even 8 days after the onset of an outburst. Their 
study showed a similar contribution to the continuum emission, as that measured during the 
peak of the outburst but also a retrograde azimuthal rotation of 58° compared to outburst 
maximum. 



Baptista et all (l2005h performed a re-ana lysis of previous spectroscopic data of IP Peg 
during outburst (IMorales-Rueda et al.l . l2000l) . based on the eclipse mapping technique. The 
difference in their new study was based on the fact that they separated the inner and outer 
regions of the emission Une light curves, obtaining in this way the "core" and "spiral arms" 
lights curves, respectively. In this way they improved the previously noisy view of the spiral 
arms and were able to trace them in detail. The two-armed spiral structures were clearly 
visible in all their eclipse maps, albeit the "blue" arm appearing farther out from the disc's 
centre than the "red" one. However, their crucial findings were the sub-Keplerian velocities 
along the spiral structures, as well as the clear correlation between the opening angle of 
the spirals and the outburst stage (the later in the outburst, the smaller the opening angle of 
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Figure 1: AAVSO long-term light curve of visual validated data. The anow shows the 
timing of our observations, just 1 day after the peak of the outburst. 



the spirals). Both these findings greatly favour the interpretation of the spiral asymmetric 
structures as tidally-induced spiral shocks instead of mere radiation patterns on an otherwise 
unperturbed AD. 



A study performed by ISteeghsl (|200lh presents a review of the spiral structure that far, 
and plots in detail the azimuthal ang les and extent of the spirals for the Hell A4686 emission 
line during an outburst presented by lHarlaftisI (119991) . In more detail, it is shown that the two 
spirals can be traced for almost 180°, while their velocity varies from 495 to 780 km s^^. 
Under the assumption of Keplerian velocities, they show that the arms cover a substantial 
part of the disc, 0.3-0.9 times the distance to the inner Langragian point Li . It is pointed out 
that tracing the spiral arms and plotting their position and intensity as a function of azimuth 
throughout an outburst, or even better applying this procedure to a range of emission lines 
simultaneously, can yield important results on the evolution of the spirals. A future paper, 
as stated in the Discussion section, is intended to address this issue. A first approach on 
tracing the spiral arms of two different emission lines (H a, Hel A6678), as well as plotting 
their emissivity as a function of azimuth, is presented in lSteeghsl (|200lb . It is shown that, 
even though there is a similar position for the spirals in terms of velocity as a function of 
azimuth, the intensity modulation is significantly different. 

Our work is part of a project of studying CVs through high resolution time-resolved 
echelle spectra and applying the indirect imaging technique of Doppler tomography. The 
simultaneous study of several emission lines can yield importa nt results on the AD structure, 
as presented for another CV system in IPapadaki et all (|2008r) . Following the introduction 
on IP Peg one can conclude on the benefits of applying such an analysis to this object, 
especially since it was observed during outburst. In this way we can study the behaviour of 
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the spiral arms and map the structure of the AD in detail. The application of Modulation 



Dopp ler tomography, an extension of the classic Doppler tomography method (ISteeghs , 



20031) . gives us the opportunity to map emission sources which can vary harmonically with 
the orbital period. Such an analysis is a first for IP Peg. 

This paper is organised as follows: in Sect. 2 we present our observations, in Sec. 3 
the data reduction, in Sect. 4 we perform the spectral analysis and Sect. 5 deals with the 
application and results of Modulation Doppler tomography. Last, Sect. 6 is devoted to the 
discussion and summary of our results. 

2 Observations 

The observing run was conducted on the 4'^ of August 1999, at ESO's NTT. Echelle spec 



troscopy was performed using the REMD mode of the EMMI spectrograph (IDekker et al. 



19861 ). The echelle grating #9 and cross-disperser #3 resulted in 21 orders covering the 
wavelength range 4000-7500 A. The dispersion increased between 0.09 and 0.17 A pix^', 
from blue to red, respectively. IP Peg was observed for 5.7 h, covering 1.5 binary orbits. In 
total, 41 spectra with an exposure time of 300 sec each, were obtained. 

Fig. [T] shows a long-term light curve of our system as generated from the AAVSO light 
curve generator tool' , choosing only visual-vahdated data to appear, and a window of a few 
days around outburst. An inspection of the outburst maximum reveals that IP Peg was at 
the peak of its outburst just 1 day before our observing run, when it was also 1 magnitude 
brighter. 

3 Data reduction 

The resulting spectra were reduced using the standard echelle reduction procedures in 
//M/l^ and the slit spectra reduction task doecslit. Whenever, both an interactive and 
non-interactive mode existed for a task, the interactive one was always chosen. 

Initially, the bias and lampflat frames were combined into one bias and flat frame with 
the zerocombine and flatcombine tasks, respectively. The spectra were then de- 
biased and a normalised flat was produced with the task apf latten in order to apply 
the flat-fielding process. Apf latten, after finding and tracing all apertures (defined as 
the sum of dispersion lines containing the signal of an order), also models both the profile 
and overall spectrum shape and removes it from the flat-field leaving only the sensitivity 
variations. The resulting normalised flat, if plotted perpendicular to the dispersion gives al- 
ternating smooth and noisy sections, smooth sections corresponding to the regions between 
the orders and having a value of 1 , and noisy sections coiTcsponding to the regions inside 
the orders and having values fluctuating around 1. The spectra were then divided by the nor- 
malised flat and thus flat-fielded. At this stage, the spectra were also corrected for cosmic 
rays with the use of the cosmicrays task, part of the crutil package. 



'|http://www.aavso.or g/datayicg/| 

^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF. 
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The task doecslit connects and combines several individual tasks to provide a single, 
complete data reduction path. It provides a degree of guidance, automation and record 
keeping. Additionally, the user with the choice of the interactive mode (used in this set of 
spectra) and the use of a graphical interface, is able to be in conti'ol of the whole procedure, 
and if necessary, fine-tune the reduction process while it runs. The several steps involved in 
doecslit and therefore appUed in the dataset are summarised as follows: 

• In an aperture reference image, preferably a bright object, the orders were located and 
given a default fixed aperture width, which could be changed interactively if needed. 

• The background sample regions were inspected. The function fitting parameters 
could be set interactively in order to best subtract the background later on. 

• Apertures were then defined for all object and standard star images, using as initial 
apertures, those defined from the aperture reference image. All images sequentially 
passed through this step. Possible deviations between the aperture reference image 
and the rest of the images, were fixed interactively by re-tracing and re-centering the 
apertures of the images. 

• The actual extraction of the spectra was done using the variance option in the weights 
parameter of the sparams task of the echelle package. Variance weighting, of- 
ten called optimal extraction, is performed by summing the pixel values across each 
aperture at each point along the dispersion. Each pixel in the sum is weighted by 
its estimated variance based on a spectrum model and the detector readout noise and 
gain parameters. For detailed inf ormation on th e optimal extraction method the reader 



should refer to Home (1986) and lMarshl (|l989f ). whose work was followed by IRAF. 



• Dispersion correction was then applied through the use of Thorium- Argon arc lamp 
spectra. Initially, a few aix lines were identified in each order (beginning, end and 
middle of the dispersion axis), making use of the Th-Ar line atlas provided by ESO 
for the specific EMMI setup used in the observing run. Additional lines were then 
automatically identified with the use of an identification Th-Ar line list available in 
IRAF. A dispersion function was then fit to the total of identified lines in each order, 
and then applied to the coiTcsponding orders of the rest of the images. The typical 
root mean square of the fits was RiO.012 A. 

• For the flux calibration we used spectra of the spectrophotometric standard LTT 7379. 
First, the bandpass regions were interactively defined and then the standard calibra- 
tion data of LTT 7379, available in IRAF, were interpolated to the user-defined band- 
passes. The flux values for each aperture were then fit with a smooth curve, in order to 
relate system sensitivity to wavelength for each aperture, thus resulting in a sensitivity 
curve for every aperture. 

• Finally the standard star spectra were extinction corrected (using the extinction file 
provided by ESO for the location of the La Silla Observatory), and flux calibrated 
using the derived sensitivity curves. 
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Table 1 : Chai'acteristics of emission lines. 
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Hel A7065 


1818 
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As a last step in doecslit, all previous steps were applied to all object spectra, re- 
sulting in the final flux-calibrated spectra of IP Peg. 

4 Spectral analysis 

The average out-of-eclipse (phases 0.15-0.85, where phase is that of mid-eclipse) spec- 
trum is shown in Fig. [2] To check that spectra of incomplete orbits do not affect the line 
profiles of the resulting average out-of-eclipse spectrum we first computed and compared 
the average out-of-eclipse spectrum of the whole dataset (1.5 orbits) and that of one full or- 
bit. The spectrum resulting from the whole dataset differed only in the B aimer lines which 
showed an enhancement of the blue-peak intensity. The rest of the lines showed minimal 
or no differences at all. Therefore, the out-of-eclipse spectra of one full orbit were used in 
Fig. 121 The orbital phases of the individual s pectra were calculated according to the linear 



orbital ephemeris given by lWolf et all (|l993l) : 



ro[HJD] =2445615.4224 + 0.'^15820616E (1) 

where Tq is the time of mid-eclipse. 

The outburst spectrum (Fig.|2]) is characterised by several double-peaked emission lines, 
along with absorption features. There is strong Balmer and high-excitation emission from 
Hell A4686 and the CIII-NIII blend. Also present is weaker emission from Hel (A4143, 
4388, 4471, 4713, 4922, 5016, 5048, 5876, 6678, 7065 and 7281), Hell (A4199, 4541 and 
5411) and CII A4267 lines. IP Peg also shows emission from the heavier elements of Fell 
(Fell 42 triplet and Fell A5169, A5303), Till A4418 and Mgll A4481. Finally, there are the 
NaD doublet and TiO absorption bands. 

All emission lines appear with symmetric line wings, except the lines of the blue part of 
the spectrum where one cannot tell with certainty due to their severe blending with neigh- 
bouring lines. Concerning the blue-to-red peak strength, the Balmer lines seem to differen- 
tiate from the high-excitation lines. They have an equal blue-to-red-peak strength (except 
Hj which has a stronger red peak) while the high excitation lines (blended or stand-alone) 
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Figure 3: Radial extension of all emission lines. Bottom x-axis gives the radial extension 
as a factor of /?wd. while top x-axis gives it as a factor of the volume radius of the primary 
Roche lobe R^ 



appeal- to have a stronger blue peak. For the Hel lines it is harder to distinguish between 
the strength of the peaks, as they are either located in the wings of stronger lines (A4143, 
A4388, A4713) or blended with equally strong lines (A4471, A4922, A5016, A4143). How- 
ever, the stand-alone Hel lines (A5876, A6678, A7065) point to a similar behaviour with 
the B aimer lines. 

In this study, we will focus on the 9 most prominent emission lines i.e. the Balmer 
lines. Hell A4686 and Hel (A 4922, 5876, 6678, 7065). Tablefflshows the characteristics of 
the emission lines (full width at zero intensity -FWZl, equivalent width -EW, and integrated 
flux). However, caution should be drawn in the interpretation of these values for the severely 
blended lines, namely H5, Hy, Hell A4686 and Hel A5016. The Balmer decrement from 
Ha to H5 1:0.77:0.68:0.55 shows a mix of optically thick and thin conditions. 

Assuming a Keplerian velocity field and AD location of line emission, we can then 
estimate the radial extension of the emission lines given the fact that the velocities measured 
from the double-peak separation and FWZl are located in the outer a nd inner disc radius. 



respectively. Using equations (2) and (3) from iPapadaki et all (120081) and taking the tidal 
effect into account the resulting radial extension of the emission lines is shown in Fig. [3] 
The bottom x-axis is measured in terms of the WD radius /?wd. while the top x-axis in 
terms of the volume radius of the primary Roche lobe 7?l- The system par ameters used to 



derive the inner and outer radius, /?in and Roau are the ones published by iBeekman et al. 



(|2000l ). The paramete rs j?wd and j?t were calculated according to equations (2.82) and 



(2.4c), respectively, in IWarnen (l2003h . A look at Fig. [3] reveals a variety of behaviours. 
The Balmer line emission commences from close to the WD, /?!„ ~ 5 ib2/?wd and extends 
up to Rout ~ 50/?wd for Ha. Both /?in and /?out gradually increase from Hj8 to Ha, albeit 
much more abruptly in Rout- The high-excitation line of Hell A4686 also appears from 
close to the WD but is constrained to a smaller outer radius /?out ~ 17/?wd and therefore 
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Figure 4: IP Peg trailed spectra of tlie 9 most prominent lines. The spectra have been 
ordered according to phase. The size of the x-axis small and big tickmarks are 35 and 70A, 
respectively. 
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Figure 5: IP Peg spectrograms. The x-axis represents the velocity in km s^' and the y-axis 
the orbital phase. The grayscale contrast has been adjusted to emphasise between emission 
(black) and absorption (white). 
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is the most centrally located line. The Hel A7065, 5876 lines have a similar /?;„ ~ 327?wd 
and 7?out ~ 35/?wd but Hel A4922 appears more centrally located, most probably due to its 
blending with the Fell line which also causes an asymmetry of the blue-to-red peak not 
present in the rest of the He lines as seen in Fig.|2] However, Hel A 6678 displays a different 
behaviour, being more centrally located and having a greater radial extension of f» 10/?wd5 
compared to the rest of the Hel Unes that are more confined with an extens ion of ^ 3/?wd- 
Non-uniform He distribution was also reported by IPiche & Szkodyi (|l989f) . when IP Peg 
was observed near outburst maximum. 

After normalising the spectra, and choosing the regions of interest, the trailed spectra of 
the 9 emission Unes were computed. They are ordered according to orbital phase and shown 
in Fig. m In order to enhance the emitting regions, the con^esponding spectrograms are 
presented in Fig. |5l with black corresponding to maximum emission. They have been folded 
on phases 0-1 and repeated over two cycles. Both figures show characteristics typical of IP 
Peg: (i) prograde rotation of a disc-like emitting region, with the blue peak being eclipsed 
prior to the red one, (ii) double-peaked emission lines signature of an AD (ii) highly phase- 
dependent double-peaked separation, (iii) an S-wave moving in anti-phase with the peaks, 
signature of chromospheric emission from the secondary star. Additional features, seen 
around ^ 1000 kms^' in the spectrograms of some emission lines, are caused by blendings 
with weaker emission lines in the red-wings of the principal lines. Each spectrogram has 
been plotted in a scale that best enhances its emitting regions, and therefore the intensities of 
the different emission lines should not be directly compai^ed. At this point, the application 
of the modulation Doppler tomography, presented in the next section, will permit us to gain 
further insight of the AD structures. 



5 Modulation Doppler tomography 



Modulation Doppler tomography, developed by lSteeghsl (l2003b . is an extension of the clas- 



sical Doppler tomography which relaxes one of its fundamental axioms, that the flux from 
any point is constant in time. This allows the mapping of time-dependent emission sources, 
and this variability information is not considered in standard Doppler tomography studies. 
Observations show that emission sources can modulate in time and this can happen either 
because there is anisotropic line-emission or because all sources are not equally visible at 
all times due to the system's geometry or both. Modulation Doppler tomography gives 
three output images; one describes the average non-varying flux distribution and the other 
two describe the sine and cosine term of the harmonically varying emission sources as a 
function of orbital phase. For a detailed description of Modulated Doppler tom ography. 



the im plementation of the MODMAP code and test reconstructions please refer to ISteeghs 



(|2003h and references therein. 

The resulting images after running MODMAP, for the 9 most prominent lines, are 
shown in Fig. |6] From top to bottom and from left to right each panel shows: the observed 
and predicted data (the latter calculated from the constant and variable components repre- 
sented by the maps), the non-variable emission distribution and the total variable emission 
component and finally the individual sine and cosine terms of the variable component. The 
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Figure 6: In each panel from top to bottom and left to right: observed data, predicted data, non-vaiiable emission distribution, total 
variable emission component, cosine and sine terms of the variable component. From left to right: H5, Hy, Hell A4686. 
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input and predicted line profiles are plotted in the same scale, reflected by the colour wedge 
on their right. The average non-varying and total modulated emission maps are plotted in 
their own scale, with the colour wedge reflecting the scale of the modulated map. The scale 
is given as a percentage of the average emission to easily depict how strong the modulation 
is, relative to the constant component. The cosine and sine amplitude maps are also plotted 
in their own (common) scale, reflected by the colour wedge on their right, again as % of the 
average emission, in order to ascertain their relative amplitudes. The Roche lobe shape and 
ballistic stream trajectory have also been computed by the MODMAP code, for the mass 
ratio and binary separation of IP Peg, and are over-plotted for reference. 

The average non-varying-emission maps of all lines show the characteristic spirals plus 
donor star emission that is well known for IP Peg. The Balmer lines show strong modulation 
of the donor star, mostly attributed to its cosine term, while Ha also shows modulated 
emission neai^ the stream trajectory. The modulation of the donor star is at a 40% level of 
the average emission in H5 and decreases when moving to higher Balmer series reaching 
about 1% for Ha. There appears to be some modulated disc emission in Ha though only 
at a 1-2% level. No modulated flux, at the signal-to-noise ratio (S/N) of our data, is seen 
from the spiral aims. 

The Hel Unes reveal a strong component at the position of the secondary star which is 
of course modulated. Two disc regions appear to show some modulation, though only at the 
1-2% level. They also show a modulated component at the inner Lagrangian point, albeit of 
low level compared to average emission. The only exception is Hel A4922 (see the cosine 
map) but one should note that this line is blended with Fell relatively strong in our set of 
spectra. Again, no significant amount of modulated flux from the spiral arms was needed to 
reproduce the observed data. 

Finally, the high-ionisation line of Hell A4686 also shows the characteristic spiral ai^ms 
and secondary star emission like the rest of the lines and similarly to the maps Doppler 



maps presented by iHarlaftisI (| 19991 ). The nearby Bowen blend causes some contamination, 
most likely the cause of the modulated structure. It is the only line in which there is a hint 
of modulated emission from the lower spiral arm at a 6% level. 



6 Summary & Discussion 

This study presents the first time-resolved echelle spectra of IP Peg, one day after the peak 
of an outburst. Our spectra show the ensemble of the characteristics previously seen in 
spectroscopic studies of IP Peg, albeit not covering all lines together. The spectroscopic 
results, deduced from the average out-of-eclipse spectra are analytically given in Sec. 5 and 
can be summarised as follows: 

• Strong Balmer double-peaked emission arising from the AD. The lines show an 
equal blue-to-red-peak strength, as opposed to the higher excitation lines that have 
a stronger blue peak, therefore pointing to a less uniform emission. The Balmer 
decrement shows a mix of optically thin and thick conditions. The radial extension of 
the emission originates close to the white dwarf and extends up to sa 0.45/?l for Ha. 
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The outer radius of the radial extension gradually increases when moving from H5 to 
Ha. 

• Strong double-peaked emission from the high excitation line of Hell A4686 with a 
stronger blue peak. It is the most centrally located emission with its outer radius 
reaching ^O.l4Ri^. 

• Weaker Hel and heavier elements (Fell, Till and Mgll) double-peaked emission. Out 
of the four Hel hues (A 7065, 6678, 5876, 4922), Hel A6678 is the most centrally 
located and shows the greater radial extension. 

The trailed spectra and spectrograms, reveal the prograde rotation, the double peaked 
emission lines arising from the disc, an S-wave component indicative of chromospheric 
emission from the secondary star and highly phase-dependent double-peaked separation, 
indicative of a distorted AD resulting in two asymmetric emission components, the spiral 
waves. All of these features, are present and obvious in the constructed Doppler maps. 

We have applied to our data set, the modulation Doppler tomography method, which 
in addition to the classic Doppler maps reveals if any harmonically varying components 
are present in the system configuration. This is of great importance in sight of the spi- 
ral waves, which nevertheless have not shown any significant harmonically varying emis- 
sion and therefore we can not claim anisotropic emission. The only line to show a hint of 
anisotropic emission, is the Hell line, and only in the lower spiral arm and at a 6% level of 
the average emission, but this could be due to the blending of the line. However, the fact that 
for the majority of lines no modulated emission is seen from the spiral arms does not mean 
that the spirals do not modulate at all, but that no significant modulated contribution was 
needed to reproduce the observed data. At higher 5/A'^ we may have seen some modulation. 
But, since we see features of amplitude a few %, this means that our S/N rules out strong 
anisotropy in disc emission. We could therefore conclude that, if present, the modulated 
emission of the spiral arms is relatively weak at that our S/N was good enough to put a 
lower detection limit of any modulated emission at 5-6%. 

The B aimer and Hel lines show strong secondary star emission, part of which is mod- 
ulated. Its level of emission decreases when moving from H5 to Ha. As stated in the 
previous section secondary star emission is naturally modulated with the cosine term. This 
is normal since irradiation induced emission from the secondary is phase depend ent, so it 



is natu ral to see a flux variation when the phase is changing. Given equation 3 of lSteeghs 



(120031) and that at phase 0.5-0.7, the face on secondary should emit more, then the cosine 
term is the dominant component. 

All Hel lines show modulated emission in 2 regions of the AD, but are at a very low 
level. Only in A4922 do they appear more pronounced. These AD modulated regions, seem 
to appear in the "gaps" between the spiral arms and modulate with the sine-term. This 
could be in accordance with the spiral arms shielding the disc areas just behind them and 
hiding a possible cosine-term of disc modulated emission as well. Moreover, the sine-term 
is expected considering their phasing is in between t he spi ral arms. Geometric shielding 
caused by the spirals was also pr oposed from Steegha (120031) and shielding of the secondary 



star by the AD was proposed by IWatson et al.l (|2003h . Finally, most of the Hel lines show 
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modulated emission near the inner Lagrangian point also covering a small region at the 
beginning of the stream trajectory. Ha though, shows pronounced modulated emission for 
the bigger part of the stream trajectory. 

In a future paper it is our aim to study the same dataset, but concerning the spiral arms 
of the different lines only. Tracing, providing the amplitude, and the contribution to the total 
flux of the spiral arms of all lines, as well as finding the ratio of the brightness between the 
two spiral arms can yield valuable results for this unique set of time-resolved spectra. 

7 The data 

• fig3.dat: an ascii file that gives the radial extension of the emission lines as a factor 
of both the white dwarf radius /?wd and the volume radius of the primary Roche lobe 
/?L. as seen in Fig. 3. 

• HJD_phase.dat: an ascii file giving the heliocentric Julian date and the orbital phase 
of the 41 spectra. 

• spectra: a folder that contains, 

- the average out-of-eclipse spectrum of Fig. 2 as a gzipped fits file, named 
lPPeg_avr.fits.gz and 

- the normalised trailed spectra of the 9 most prominent emission lines as shown 
in Fig. 4. These trailed spectra were the ones used for the reproduction of Fig. 4- 
Fig. 6. The naming convention is {emission line}_{sequence number}. fits.gz, 
resulting in 9 (number of emission lines) x 41 (number of spectra) different files. 
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